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^ ' We report the discovery of a single-lined spectroscopic binary with an Ultra 

^ ' Cool Dwarf (UCD) primary with a spectral type between M8 and L0.5. This sys- 

^ ■ tern was discovered during the course of an ongoing survey to monitor L dwarfs for 
radial velocity variations and is the first known small separation (a < 1 AU) spec- 

^ ■ troscopic binary among dwarfs at the M/L transition. Based on radial-velocity 

■ measurements with a typical precision of 300 m s~^, we estimate the orbital 
^ ; parameters of this system to be P = 246.73 ± 0.49 d, ai sinz = 0.159 ± 0.003 

I ■ AU, M2 sin i = 0.2062(Mi + M2)2/3 ± 0.0034 Mq. Assuming a primary mass of 

'sj" ■ Ml = O.O8M0 (based on spectral type), we estimate the secondary minimum 

^ ' mass to be M2 smi = O.O54M0. With future photometric, spectroscopic, and in- 

O ■ terferometric observations it may be possible to determine the dynamical masses 

■ of both components directly, making this system one of the best characterized 
^ ' UCD binaries known. 
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Introduction 



Near infrared (NIR) surveys, such as 2MASS, DENIS, and SDSS (witli its 2;-band ca- 
pability), have resulted in a rapid increase in our knowledge of the properties of stars. This 
is particularly t rue for the late M , L and T spectral types, collectively known as Utracool 
Dwarfs (UCDs, Kirkpatrick 2005 ). Today we know of more than 600 L and T dwarf^. De- 
spite these discoveries, we still do not have a clear understanding of how molecular cloud 
material assembles itself into such relatively low mass objects and fundamental properties 
such as temperature, age and mass remain coarsely determined. Binary star systems are 
a crucial tool for addressing both of these issues. UCD binaries allow us to directly mea- 
sure the masses, and possibly the radii, of these objects and constrain theoretical models 
of the structure and emergent flux (i.e. IStassun et al.ll2006l ). Even without determining 
absolute ages, assuming that the components of a binary have identical ages can provide 
constraints on atmospheric models through estimates of the luminosity ratios. Different 
formation scenarios make varied p redictions for the statistical properties of UCD binary sys- 
tems (see iBurgasser et al.ll2007al for an overview). As a result, studying the properties of 
UCD binaries may shed light on the formation mechanism of the entire class of objects. 

Although models of field L and T dwarfs have reached a high degree of sophistication, 
model spectra, particularly for L dwarfs, are highly dependent upon the assumed cloud 
model. ICushing et al.l (120081 ) . for example, found that changes in the assumed degree of cloud 
sedimentation could alter the derived effective temperature, Tes, by up to several hundred 
Kelvin and the log of the surface gravity, logg, by 0.5 dex. A spectroscopic binary offers 
the prospect of constrained masses and coevality and would provide exce llent simultaneous 
const raints on the cloud model, object masses, and effective temperatures (IMarley fc Leggett 
2008h . 



The binarity of UCDs has been studied both spectroscopica lly and with direct imaging. 
Relative to main sequence stars (i.e. iDuquennoy fc Mayoiill99ll ). the UCD binaries are more 
rare, lie at closer separ ations, and are more likely to have components with equal masses 
(IBurgasser et al.ll2007al and references therein). Owing to the difficulties in detecting close, 
high contrast ratio binaries, imagin g surveys do not provide a clear picture of U CD binarity 



at sep a rations less than 1 AU (see iReid et al 



20081) . D espite these difficulties, iBouy et al. 



(I2OO4J ). iLane et al.l (1200 ll ). and iGolimowski et al.l (120071 ) have used high-resolution imaging 
to determine the orbit of young L or M dwarf binaries. The binaries with small separations 
are more readily detectable by searching for the radial velocity signal due to the reflex motion 
of the primary star. Spectroscopic searches for UCD binaries have identified several binaries 
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and candidate binaries (|Joergens fc Miillerl 120071 : iKenyon et al.l l2005l : iGuenther fc Wuchterl 
20031 : iBasri fc Reinersll2006[). though t he tota l number of observations of each object tends to 
be small. Recent work by iBurgasserl (l2007bl ) has demonstrated that low-resolution spectra 
alone can be used to detect the faint companions by searching for the subtle signature of the 
cooler object at the wavelengths of certain molecular features. 

We presented results fr om a pilot study ta rgeting nine field L dwarfs with the Phoenix 
instrument on Gemini-S in iBlake et al.l (120071 ). Here, we report the detection of a single- 
lined UCD spectroscopic binary (SBl) with an orbital separation of approximately 0.4 AU. 
2MASSJ03202839-0446 358 (hereaf t er 2M 0320-04) was included in the catalog of nearby, 
cool stars presented by I Cruz et al.l (120031) who spe c trosco pically classified it as a possible 
M8 dwarf based on optical spectra. I Wilson et al.l (120031 ) classified the object as an L0.5 
dwarf based on low resolution NIR spectroscopy. In §2 we describe the high-resolution NIR 
spectroscopy of 2M0320— 04, in §3 we describe the modeling process used to extract radial 
velocities from these data, the fit of a Keplerian orbit to the radial velocities, and the search 
for the spectral lines of the secondary, and in §4 we describe the implications of this discovery 
for future studies of UCDs. 



2. Observations 



We observed 2M0320-04 with the NIRSPEC spectrograph (IMcLean et al.l [19981) on 
the Keck telescope on 14 nights between September 2003 and January 2007 as part of a 
program to monitor a large sample of L dwarfs for radial velocity variations. In all cases, 
we used the same instrumental setup. This setup was selected to cover the CO bandhead 
and R-branch features around 2.3/im. This spectral region is rich in telluric absorption 
features due to methane. As described in the next section, these telluric features will serve 
as the wavelength reference for our radial velocity measurements. We used a 3 pixel (0.432") 
slit, with the N7 blocking filter, to produce an approximate scale of 0.3 A pixel and a 
resolution of R = A/AA ~ 25000. The extracted spectra contain 1000 pixels and cover the 
approximate spectral range 2.285/im to 2.318/im. Exposures ranged from 900s to 1200s per 
nod position with adjustments made according to observing conditions in order to maintain 
approximately constant S/N. The data were gathered in nod pairs so as to facilitate the 
subtraction of sky emission lines. In total, we collected 16 nod pairs. 
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3. Analysis 



After subtraction of the nod pairs to remove sky emiss ion line s , we e xtracted the spec- 
tra following the optimal extraction procedures outlined in iHornd (119861) . We r aodel ed the 
extracted spectra following a procedure similar to that described in lBlake et al.l (120071 ). Our 
model begins with two high-resolution template spectra: one for the Earth's atmosphere and 
one for the L dwarf. The high-resolution (5 x 10~^u m pixel"^) spectrum of the Earth's atmo- 



sphere is provided by [Livingston fc Wa 



lacd (119911). The hi gh-resolution synthetic spectra of 



L dwarfs were computed as described in lMarley et al.l (120021 ). with a number of improvements 
to be described in a fu t ure p ublication. The models apply the condensation cloud model 
of lAckerman fc Marleyl (120011 ) with a sedimentation parameter of fsed = 3, corresponding 
to a moderate amount of condensate settling. The models used her e have solar metallicity 
( lLoddersll2003l ). use the opacities described in lFreedman et al.l ( 120081 ) . and a fixed gravity of 
logg = 5 (cgs) and cover a range of Tes from 1200 to 2400 K. The synthetic spectra provide 
monochromatic fluxes spaced 4.2 x 10~^/im apart. We convolve and re-sample the product 
of the telluric and L dwarf high-resolution spectra to generate the model that we then fit to 
the extracted 1-D spectra. Our model has several free parameters. The parameters related 
to the L dwarf are the projected rotation velocity of the L dwarf (V^sinJ) where / indi- 
cates the inclination of the rotation axis to the line of sight, its T^q, and its radial velocity. 
The parameters related to the spectrograph are the PSF width, flux normalization, and the 
wavelength solution (i.e. the mapping from wavelength to pixel position). 

For the parameters of the L dwarf we first determined a best initial fit with V sini = 
16.5 ± 0.5 km s~^ using the L dwarf spectral template with T^g = 2200K using least squares 
fitting. This spectral tem plate provided the best overall fit to the da ta and is also consistent 
with the LO spectral type (IBasri et al.ll2000l : iGolimowski et al.ll2004l ). We modeled the wave- 
length solution as a third-order polynomial, the overall flux normalization as a third-order 
polynomial, and the spectrograph PSF as a single Gaussian. The PSF is expected to be 
asymmetric at some level, and to vary across the spectrograph order. Spectral models that 
accommodate these subtleties may yield even higher precision than that achieved here. With 
the Vsinl and T^s fixed the model has a total of ten free parameters. For each spectrum, 
we determined the best fit values for each parameter, including the radial velocity, using the 
AMOEBA algorithm. An example of a NIRSPEC spectrum and best fit model is shown in 
Figure 1. The algorithm failed to converge on a solution for one out of 28 spectra, resulting 
in a total of 27 radial-velocity measurements. Visual inspection of this spectrum indicated 
a strong flux gradient across the order, likely due to a poor extraction. We used a bootstrap 
simulation t o estimate the erro r on the individual radial velocities, similar to the technique 
described in lBlake et al.l (120071 ). We simulated NIRSPEC observations of L dwarfs with dif- 
ferent S/N and radial velocities and used our modeling procedure to estimate our ability to 
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recover the known radial velocity. Since the NIRSPEC data are of relatively high S/N 60 
per pixel), and the VsinJ is a modest 16.5 km s~^, the results of our simulation indicate 
that the statistical errors on the radial velocities should be rather small (~ 60 m s~^). The 
measured radial velocities are listed in Table 1. 

An initial search for periodicity in the radial velocities with the Lomb-Scargle Peri- 
odogram revealed a clear signal with a period P ~ 250d, prompting a more detailed anal- 
ysis. The radial velocities covering nearly 5 cycles of the binary were fit with a Keplerian 
orbital model using standard non-linear least-squares techniques. The six parameters of this 
model are the period, P, the systemic velocity, 7, the radial velocity semi-amplitude, Ki, 
the eccentricity, e, the time of periastron passage, Tq, and the longitude of periastron, u. 
The scatter in the fit was found to be significantly larger than the internal uncertainties, 
possibly as a result of systematic errors. For the final fit we adjusted the internal errors by 
adding 0.32 kms~^ in quadrature to the error estimates based on the bootstrap simulations, 
so as to force a reduced value near unity. The resulting orbital elements are listed in 
Table 2, wherein we also state our estimate of the projected semi-major axis of the orbit 
of the primary, aisini, and the observations and orbital fit are shown in Figure 2. The 
eccentricity (e = 0.065 ±0.016) is small but significant at the 4a confidence level. The radial 
velocity semi-amplitude of the primary, Ki = 7.02 ± 0.12 km s~^, is 21 times the per-point 
measurement precision. 



4. Discussion and Conclusions 



Based on the estimated spectral type (M8 to L0.5) we can estimate the effective tem- 
perature o f the primary of 2M0320 — 04 if we assume that the light from the secondary is 
neghgible. iGolimowski et al.l (120041 ) present estimates of T^s as a function of spectral type 
and from this work we estimate 2200 < Tcs < 2400K for an assumed age of 3 Gyr, in gen- 
eral a greement with our fits to the L dwarf templates. Using the models of iBaraffe et al. 



( 120031 ) to relate T^s to mass at an age of 3 Gyr, we estimate that the primary has a mass 
0.075 < Ml < O.O8IM0. We can use the mass function derived from the radial-velocity 
fitting procedure to estimate the mass of the secondary as a function of sini where i is the 
inclination of the system to our line of sight. If we assume Mi = O.OSM fr^, and that sini ^ 1 , 
then the mass of the secondary is Mp, = .054Mf7i. Using the models of iBaraffe et al.l (120031 ) 



and the work of IGolimowski et al.l (120041 ). the secondary would have Tgfj ~ 1350K with a 
spectral type of approximately L7 to T3 at the assumed age of 3 Gyr. As sini decreases, the 
total mass of the system increases, resulting in an observed spectrum that is not dominated 
by the LO primary. This would be inconsistent with observations without requiring that the 
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more massive component of this system be significantly under- luminous. In particular, val- 
ues of sini < 0.7 would result in a secondary star that is more massive, though less luminous, 
than the primary. While it is possible that sinz ~ 1 and the system presents eclipses, the 
relatively small value of -|- R2)/cl ~ 10~^ makes eclipses unlikely. Given the constraint 
sini > 0.7 the actual probability of observable eclipses is likely somewhat higher. 

While we have no direct measure of the age of the system, it is possi ble to constrain th e 
age from its kinematics. Using the distance es timate of 26. 2± 4.3 pc from I Cruz et al.l (120031 ). 
the proper motion of 0.678 ± 0.038" year~^ (jPeacon et al.ll2005l ). and our measurement of 
the radia l velocity, we can calculate the space velocity and U,V,W velocity components 
following I Johnson fc Soderblonj (119871 ). Since the parallax of this system is not known, the 
large error on the distance estimate results in large errors on the kinematic estimates. We 
find velocity components [U,V,W] = [-62.0±9.3,-38.6±8.9, -34.3±7.2] km s"^ r elativ e to 
the local standard of rest. Using the age- velocity relation from Eqn. 8 of IWieleru (119771 ) we 
estimate an age based on the W velocity of r > 3.6 Gyr (95% confidence). 

It is interesting to consider the detectability of the spectral lines of the secondary in our 
spectra, whic h woul d allow us to estimate directly the mass ratio. Based on the models of 
Baraffe et al.l (120031 ) . the expected i^-band flux ratio for objects of these masses is ~ 10 at 
T = IGyr and ^ 100 at r = 5Gyr. While detecting the secondary lines at the later age would 
be challenging, methods like TODCOR (IZucker &: Mazehlll994l ) have been used to recover 
secondaries in systems with flux ratios ~ 50 (D. Latham; private communication). Following 
a method similar to TODCOR we searched for the spectral lines of the secondary as follows. 
Using the orbital solution presented in Table 2, we searched a grid in two parameters; the 
mass of the secondary and the i^-band flux ratio. At each grid point, a second theoretical 
template spectrum with T^fi = 1400K was added to the fitting procedure described in Section 
3. The V sin I of the secondary was assumed to be the same as that of the primary. During 
this process the radial velocity of the primary of 2M0320— 04 was fixed to the value from 
the orbital solution, the radial velocity of the secondary was also fixed based on the orbital 
solution and the assumed value of M2, and the i^'-band fiux ratio was fixed. At each grid 
point this modified modeling scheme was used to determine the best fit of this model to the 
data. Significant improvements in from the case of a fiux ratio of 0.0 would indicate the 
detection of the secondary component. We carried out this procedure for a subset of 12 of 
our spectra gathered near times of quadrature, when the velocity separation of the primary 
and secondary would hkely be greatest. We found no evidence for the spectral signature 
of a secondary component with a fiux ratio greater than 0.1 in band. We note that the 
non-detection of the secondary spectral lines also implies sini ^ 1.0. If sini = 0.8 the 
expected fiux ratio at r = 5 Gyr would be ~ 0.3 and the secondary would likely have been 
detected. Detection of the spectral lines of the secondary and resolution of the system with 
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direct imaging would allow for the first direct measurement of the mass of a field T dwarf. 
Since the orbital solution combined with direct imaging provides a distance measurement, 
future observations of this system could also provide important empirical tests of theoretical 
models for old objects at such low masses. 

This work demonstrates the importance of radial velocity searches for binary UCDs 
with small (a < 1 AU) separations as a complement to the direct imaging searches for 
companions at larger separations. A more detailed analysis of our data, including a more 
sophisticated model of the individual spectra, may lead to the detection of the spectral 
lines of the secondary and an estimate of the mass ratio. If the inclination can also be 
measured, then the masses o f both components may be determi ned. If we assum e the 
models of iBaraffe et al.l (120031 ) . sini = 1.0, and a distance of 26 pc (ICruz et al.ll2003l ). then 
the maximum angular separation of the pair is ~ 16 mas, below the capabilities of the Keck 
Laser Guide Star AO system. While the fiux ratio may be large {AK ^ 5 mag, depending on 
age), future interferometric systems may be able to resolve both components and provide the 
measurement of the inclination required to directly measure the masses of both components. 



Note: During the completion of this Letter we became aware of work by lBurgasser et al. 



(120081 ) describing a tentative det ection of th e secon dary in the 2M0320— 04 system using the 
spectral diagnostics described in iBurgasserl (l2007bl ). 
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Table 1. RV Data for 2M0320- 



HJD-2400000 


RV 






km s ^ 


km s ^ 


52921.0960 


6.95 


0.33 


52921.1102 


6.26 


0.33 


52922.1054 


6.70 


0.33 


52922.1196 


6.35 


0.33 


52957.0218 


4.18 


0.34 


52957.0361 


4.49 


0.33 


53272.1197 


-6.17 


0.33 


53272.1304 


-6.85 


0.33 


53273.0793 


-6.56 


0.33 


53273.0900 


-5.97 


0.33 


53328.8225 


-3.97 


0.34 


53328.8343 


-4.03 


0.33 


53421.7131 


6.74 


0.33 


53421.7239 


6.18 


0.32 


53669.8797 


6.53 


0.32 


53669.8919 


5.84 


0.33 


53670.8719 


6.37 


0.33 


53670.8841 


6.26 


0.33 


53686.8538 


5.59 


0.33 


53686.8659 


5.78 


0.33 


53742.8030 


-3.25 


0.33 


53742.8151 


-2.92 


0.33 


54023.9636 


-7.25 


0.33 


54023.9757 


-7.70 


0.33 


54100.7404 


2.03 


0.33 


54100.7526 


1.03 


0..33 


5'1101.7i()8 


l.()7 


l).:-!3 



Note. — Individual radial-velocity 
measurements. Errors include the 
estimated statistical errors added in 
quadrature with a 320 m s~^ system- 
atic error. 
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Table 2. Orbital and System Parameters 



Parameter 


Value 


Units 


P 


246.73±0.49 


days 


7 


-0.063 ± 0.078 


km s~^ 




7.02 ±0.12 


km s~^ 


e 


0.065 ± 0.016 




U! 


177 ± 17 


o 


To 


2453537 ± 11 


HJD 


ai sin i 


23.75 ±0.41 


10" m 


M2 sin i 


0.2062(Mi + M2)(2/3) ± 0.0034 


Mq 


Vsini 


16.5±0.5 


km s~^ 


J 


13.259±0.024 




H 


12.535±0.023 




Ks 


12.134±0.026 




RA 


03:20:28.39 


hh:mm:ss (J2000) 


DEC 


-04:46:36.4 


dd:mm:ss (J2000) 



Note. — Derived and observed parameters of the 2M0320— 04 sys- 
tem. 
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30001 




2.29 2.30 2.31 2.32 

Fig. 1.— Example of a NIRSPEC spectrum of 2M0320-04 (bottom) along with the best fit 
model offset by a constant value (top). The spectra consist of rotationally broadened CO 
and H2O features from the L dwarf along with narrow telluric CH4 lines. 
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10^ ^ 




53000 53200 53400 53600 53800 54000 54200 

HJD-2400000 



Fig. 2. — Radial-velocity measurements and orbital solution for 2M0320— 04 plotted as a 
function of time. Our observations span approximately five orbital periods. 



